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Abstract 

The uncertainty in the absohitc value of the air-fluorescence yield 
still pTits a severe limit on the accuracy in the primary energy of ultra- 
high-energy cosmic rays. The precise measurement of this parameter 
in laboratory is in turn conditioned by a careful evaluation of the en- 
ergy deposited in the experimental collision chamber. In this work 
we discuss on the calculation of the energy deposition and its accu- 
racy. Results from an upgraded Monte Carlo algorithm that we have 
developed are compared with those obtained using Geant4, showing 
excellent agreement. These updated calculations of energy deposition 
are used to apply some corrections to the available measurements of 
the absolute fluorescence yield, allowing us to obtain a reliable world 
average of this important parameter. 



1 Introduction 



The fluorescence technique for detection of extensive air showers has been 
proved to be very fruitful. It is based on the measurement of the faint fluores- 
cence radiation of molecular nitrogen in the ~ 300 — 400 nm spectral range 
induced by charged particles (mostly electrons) of a shower. By means of 
imaging telescopes, the longitudinal profile of the shower is recorded allowing 
its geometrical reconstruction. In addition, it provides a nearly calorimetric 
measure of the energy of the primary particle by means of the air-fluorescence 
yield Y, that is, the number of photons per unit deposited energy. This 
key parameter is a function of wavelength and depends on the atmospheric 
parameters, i.e., pressure, temperature and air composition. Nonetheless 
the energy scale of fluorescence telescopes is mainly determined by its abso- 
lute value, usually given at standard atmospheric conditions (e.g., dry air at 
1013 hPa and 293 K) and for the most intense band of the spectrum, located 
at 337 nm. 

The air-fluorescence yield is measured in dedicated experiments in labo- 
ratory using electron beams to excite air at known conditions in a collision 
chamber. The induced fluorescence is registered by an appropriate optical 
system. Several experiments use electrons from a ^"Sr-^^Y radioactive source 
(average energy around 1 MeV), whereas other experiments use beams from 
accelerators or electron guns (energies ranging from keV to GeV). It has been 
shown that the fluorescence yield is nearly independent of electron energy [1}- 
|3], which is the underlying principle of a calorimetric determination of the 
primary energy of extensive air showers. In fact, it is also expected to be 
independent of the type of incident particle, since both fluorescence emis- 
sion and energy deposition are governed by low-energy secondary electrons 
produced in ionization processes. 

An accurate experimental determination of the absolute air-fluorescence 
yield requires an end-to-end calibration of the optical system as well as a 
careful Monte Carlo (MC) evaluation of the energy deposited in the field of 
view. We discuss on the calculation of the energy deposition and its impact 
on the determination of the absolute air-fluorescence yield in dedicated ex- 
periments. Results of energy deposition from an upgraded MC algorithm 
that we have developed are compared with those obtained using Geant4. 
Finally, corrections previously proposed in |lHn] to some fluorescence-yield 
measurements have been updated according to these calculations, and a re- 
liable world average of the absolute air-fluorescence yield is presented. 
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Figure 1: Comparison of results of electron energy deposition from our up- 
graded MC simulation and Geant4 for a cube 10 cm side filled with air at 
1.29 mg/cm'^. An excellent agreement (< 1%) is observed in the wide energy 
range from 1 MeV to 100 GeV. The energy loss given by the Bethe-Bloch 
formula is also shown to illustrate the contribution of escaping 6 rays. 

2 Calculation of energy deposition 

Collision stopping power is accurately described by the Bethe-Bloch theory, 
reference data being tabulated in [7] for a number of materials including air. 
However, for a correct determination of the energy deposited in a finite air 
volume, the fraction of energy transferred to secondary radiation escaping 
this volume has to be evaluated. In particular, high-energy secondary elec- 
trons {S rays) have large ranges. For instance, 100 keV electrons can travel 
14 cm in air at atmospheric pressure before being stopped. Although 6 rays 
are barely produced, they carry a significant fraction of the energy lost by 
the primary particle. This effect causes the energy deposition to be lower 
than the total energy loss by ~ 8% for incident electrons of 1 MeV and by 
more than 30% for GeV electrons at usual experimental conditions (i.e., at- 
mospheric pressure and a cm-sized collision chamber). This is illustrated in 
figure [T| where MC calculations of the energy deposition E^ep in a 10 cm cube 
as a function of electron energy are plotted together with the Bethe-Bloch 
energy loss. 

In addition, other features related to the propagation of electrons through 
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the experimental collision chamber have minor effects on the actual energy 
deposition. Firstly, stopping power varies along the trajectories of electrons 
as they lose energy gradually. In the second place, scattering due to both 
elastic and inelastic collisions affects the average path length of electrons 
inside the considered volume, contributing to either a rise or a reduction of 
the energy deposition depending on geometry. 

The energy deposition in the collision chamber of fluorescence-yield ex- 
periments is usually determined by available MC codes like Geant4 [8] and 
EGS4 [9] . We have developed an alternative MC algorithm [H [10] using our 
own compilation of cross sections and other molecular data. Unlike standard 
MC codes, which normally employ the multiple scattering approximation, our 
program simulates all the individual interactions of electrons. Comparison 
of these independent simulations is thus a valuable test for the correctness 
and consistency of results on energy deposition, allowing estimation of un- 
certainties. This is mandatory to achieve high precision level (< 5%) in the 
absolute fluorescence yield. 

3 Comparison of Monte Carlo simulations 

A series of comparisons between our MC algorithm and Geant4 has been 
reported in previous works [H-E], showing general agreement. However, non- 
negligible discrepancies were found and they have been being investigated. 
As a result of these studies, several improvements have been made in our 
MC algorithm: 

1. Corrections for the density effect in molecular cross sections, only no- 
ticeable at very high energy, have been reviewed and now are applied 
to the cross section for K-shell ionization as well (see for details). 

2. The energy distribution of secondary electrons has been modified to 
better reproduce the end tail of S rays given by Moller scattering. 

3. Angular-differential cross sections (both elastic and inelastic) have been 
corrected for some relativistic effects. 

Updated results from our MC algorithm are presented here and compared 
with results from Geant4. In figure [T| a comparison for the average energy 
deposition per electron in a cube of 10 cm side filled with air at 1.29 g/cm^ 



3 



(1013 hPa and 273 K) is shown. The geometry was defined such that elec- 
trons crossed the cube along an axis going through the center of opposite 
faces. For this comparison, the standard physics list G4emStandard_opt3 of 
Geant4 version 9.4, including multiple scattering, ionization, bremsstrahlung 
and emission of both X rays and Auger electrons, was used. An excellent 
agreement (< 1%) holds between both simulations in the wide energy range 
from 1 MeV to 100 GeV. Different geometries and other physics lists of 
Geant4 have been tested, always showing excellent agreement. 

In [6], we reported some minor disagreement (~ 2%) between our MC 
algorithm and Geant4. They are mainly attributed to the inaccuracies in the 
energy distribution of secondaries in earlier versions of our code (improve- 
ment 2). The corrections in the angular-differential cross sections (improve- 
ment 3) turned out to be appreciable only at low energy (< 1 MeV), since 
high-energy electrons follow almost rectilinear trajectories at those condi- 
tions. The latter corrections also contribute to a better agreement between 
both simulations. 

On the other hand, we still find some differences between our MC algo- 
rithm and Geant4 in the simulated trajectories of electrons, which may be rel- 
evant at keV energies or in experiments with strong dependencies on geomet- 
rical factors. In figure [2| histograms of track lengths of 1 MeV electrons cross- 
ing a 10 cm diameter sphere filled with air at atmospheric pressure are shown 
for data from our upgraded MC algorithm and from Geant4 with the stan- 
dard physics list or the low-energy PENELOPE extension (G4emPenelope 
library). Notice that data from our MC algorithm are more spread and its 
mean value is slightly lower than the predictions from both Geant4 simu- 
lations. These discrepancies could be due to the approximation of multiple 
scattering employed by Geant4 and are under study. 

Comparisons between our MC algorithm and Geant4 for the pressure 
dependence of the energy deposition have also been carried out. For this 
purpose, the energy deposited per unit mass thickness (MeV g~^ cm^) is 
calculated as the ratio of the average energy deposition and the average 
mass thickness X traveled by incident electrons. Note that the trivial linear 
growth of the total energy loss with mass thickness vanishes in this ratio 
and it solely remains the pressure dependence due to the escaping energy 
carried by secondary radiation. This parameter is also insensitive to possible 
slight differences between the simulated track lengths of incident electrons. 
In figure [3| results of the E^f.p/X ratio for 14 MeV electrons crossing a 10 cm 
sphere are shown for the 0.1 — 10000 hPa air pressure range (293 K). Very 
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Figure 2: Track lengths of 1 MeV electrons going through a 10 cm sphere 
(air at atmospheric pressure). From left to right, histograms are obtained 
using our MC algorithm, Geant4 with the standard physics list and Geant4 
with the low-energy PENELOPE extension. 

good agreement between both simulations is observed again, obtaining an 
approximated power law with a marked change of slope at around 10 hPa, 
which is associated to the average range of 410 eV photons emitted by ni- 
trogen after K-shell ionization (see |T0] for a thorough interpretation of sim- 
ulation results). In this comparison, the low-energy PENELOPE extension 
of Geant4 has been used instead, as it allows to reduce the cut in energy for 
production of secondary radiation down to a value as low as 10 eV (the lower 
the pressure, the lower the energy of secondary particles that can escape the 
sphere). The standard implementation of Geant4 gives very similar results 
at pressures above 15 hPa. Best fits of data in the 15 — 1000 hPa interval 
(typical relevant pressures for cosmic-ray detection) are also shown in the 
figure. 

It is worth mentioning that our Geant4 simulations are fully consistent 
with those performed by the MACFLY [12] and AIRFLY [llj experiments, 
which used this MC code to evaluate the energy deposition in their collision 
chambers. Some differences with respect to Geant4 calculations made by the 
AirLight experiment [11] are found, nevertheless the authors of this experi- 
ment reported very recently [15] a bug in their simulation that likely justifies 
these discrepancies. Dedicated comparisons with other MC codes have not 
been performed yet, however tests have also been made [H [5] using EGS4 
data given by the FLASH experiment for 28.5 GeV electrons [13]. These 
calculations deviate more significantly from ours (up to 8%, depending on 
pressure), although the discrepancies could fairly be explained by a possible 
inaccurate treatment in [13j of the density effect, which is very significant 
at such an energy. Geometrical factors related to the spatial distribution 
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Figure 3: Comparison between our MC algorithm and Geant4 for the pressure 
dependence of the energy deposition per unit mass thickness of air for 14 MeV 
electrons crossing a 10 cm diameter sphere. Both simulations predict a power 
law with a marked change of slope at around 10 hPa. Best fits of data in the 
15 — 1000 hPa range are also shown. 

of energy deposition in this experiment have also been investigated, finding 
some discrepancies between the simulations that are still unclear. 

Taking into consideration that our MC algorithm and Geant4 use differ- 
ent approaches for particle transport and are based on nearly independent 
molecular data, this consistency of results suggests that a high-precision level 
in the energy deposition is reached. On the other hand, this precision is lim- 
ited by the accuracy in reference data of the Bethe-Bloch stopping power, 
which is estimated to be about 1% in the energy range of interest [7j. In 
view of that, we estimate a conservative total uncertainty of 2% in the above 
calculations of energy deposition. 

4 Impact on the fluorescence yield 

Several absolute measurements of the air-fluorescence yield are available in 
the literature. In recent experiments [ISHIHIIS], energy deposition was care- 
fully evaluated by means of detailed Geant4 or EGS4 simulations. However, 
the authors of other well-known experiments piTHTU] determined it from the 
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Bethe-Bloch formula ignoring the effect of the escaping secondary radiation, 
which can cause large systematic errors, as discussed above. 

In previous works [1H6] , we proposed corrections to the energy deposition 
calculated in some of these experiments and the fluorescence-yield values were 
re-normalized correspondingly (the experimental fluorescence intensities were 
unchanged). As a consequence of these corrections, a very consistent sample 
of fluorescence-yield values was obtained, allowing us to calculate a reliable 
world average of this important parameter. The above-mentioned upgrades 
in our MC algorithm slightly affect this average. Next, updated results from 
this analysis are presented. 

Figure |4] shows a comparison of the available experimental data of the 
absolute fluorescence yield normalized to common conditions, i.e., 337 nm 
band, dry air at 1013 hPa and 293 K|^ Original fluorescence yields (after be- 
ing normalized to these reference conditions) are plotted on the left. In the 
right plot, results from experiments that did not evaluate the energy deposi- 
tion by simulation are corrected according to our calculations. Experimental 
uncertainties are not modified in any case, and it has been assumed that 
different measurements of a same experiment (only mean values are shown 
in|4]) are fully correlated but no correlation exists between experiments. The 
sample includes the very precise measurement recently reported by the AIR- 
FLY Collaboration [TB], that use Geant4 to evaluate the energy deposition 
assuming an error contribution of 2% to the total uncertainty of 4%. 

Our corrections increase the fluorescence yields of [T7I4T9] by 9 — 24%, 
leading to an improved compatibility of results. Notice that the statistic 
normalized by the number of degrees of freedom is lowered from 1.81 to 0.87 
after our corrections. The sample average (F) = 5.58 ph/MeV is calculated 
using as weights the reciprocals of the squared experimental uncertainties. 
Somewhat different averages ranging from 5.46 ph/MeV to 5.68 ph/MeV are 
derived from plausible variations on the statistical analysis (e.g., treatment 
of outliers, weights and correlations) and on our corrections (e.g., geomet- 
rical implementation details in our simulations). The standard error of the 
weighted mean (T(y) = 0.15 ph/MeV (2.6%) is determined from the quoted 
uncertainties following the usual procedure. However, this low uncertainty 
seems to us to be unrealistic, because some experiments did not include any 
error contribution from the evaluation of the energy deposition, which we 

^Quenching data as well as intensities of the remaining bands relative to the 337 nm 
band necessary for this normalization have been taken from [llj . 



7 



<Y> = 5.31 ph/MeV 
o-<y> = 0.19ph/MeV 
X^/ndf = 1 .81 



y-L (ph/MeV) 



AIRFLY 

AirLight 

FLASH 

MACFLY 

Lefeuvre 

Nagano 

Kakimoto 



<Y> = 5.58 ph/MeV 
o-<,> = 0.15ph/MeV 
Z^/ndf =0.87 



Tfi. prediction 



Corrected 



y337 (ph/MeV) 



Figure 4: Effect of our corrections on the available fluorescence-yield mea- 
surements normalized to common conditions (1013 hPa, 293 K and 337 nm 
band). Original results are on the left and corrected ones (when applicable) 
are on the right, showing better compatibility. Vertical lines represent the 
position of the sample weighted average (Y) and its standard error a^y) ■ The 
theoretical fluorescence-yield value predicted by our MC simulation is also 
shown for comparison. 

estimate to be 2% at least for Geant4 and our MC algorithm. In addition, 
there are still some discrepancies between simulations to be elucidated, as 
pointed out above. Taking into account all these considerations, we propose 
an average fiuorescence-yield value of 5.58 ph/MeV and a conservative uncer- 
tainty estimate of 4%, i.e., as low as the smallest uncertainty of the sample 
corresponding to the AIRFLY measurement. Indeed, it can be understood 
from the present analysis that this average and the AIRFLY measurement 
(5.61 ph/MeV) are complementary results showing high reliability. 

Our MC algorithm is able to simulate both the energy deposition and 
the fluorescence emission [TU], which in combination with experimental 
quenching data from [TT] also allows us to give a prediction of the absolute 
fluorescence yield of 6.3 ± 1.3 ph/MeV, in agreement with experimental data 
(figure |4| right plot). Although we estimate a large uncertainty in this value 
due to the uncertainties in the many molecular parameters involved in the 
calculation, the agreement with the experimental results provides a valuable 
theoretical support to these measurements. 
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5 Conclusions 



We show that effects of secondary radiation escaping the colhsion chamber 
are very relevant for a correct determination of the energy deposition in 
experiments dedicated to measure the air-fluorescence yield. Assuming the 
Bethe-Bloch stopping power instead of computing the energy deposition by 
simulation may lead to large systematic errors in the fluorescence yield. 

New results of energy deposition from an upgraded MC algorithm that 
we have developed are presented here, showing an excellent agreement with 
Geant4 for electrons with energies in the wide 1 MeV — 100 GeV range. 
Prom the present comparison, we conclude that these calculations have a 
precision level of 2%. On the other hand, small differences in the simulated 
electron trajectories are found, with possible implications at lower energy or 
in experiments with strong dependencies on geometrical factors. 

Our updated calculations of energy deposition have been used to correct 
experimental fluorescence-yield values reported by experiments that did not 
evaluate the energy deposition by simulation. These corrections improve 
signiflcantly the compatibility of the available measurements of the absolute 
air-fluorescence yield. As a result of this analysis, we present a reliable 
average of 5.58 ph/MeV with an estimated uncertainty of 4%. In addition, 
experimental data are shown to be consistent with a theoretical prediction 
of the absolute fluorescence yield that we obtained by means of our MC 
simulation. 
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